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We examine characteristic properties of deterministic and stochastic diffusion in low-dimensional 
chaotic dynamical systems. As an example, we consider a periodic array of scatterers defined by 
a simple chaotic map on the line. Adding different types of time-dependent noise to this model 
we compute the diffusion coefficient from simulations. We find that there is a crossover from de- 
terministic to stochastic diffusion under variation of the perturbation strength related to different 
asymptotic laws for the diffusion coefficient. Typical signatures of this scenario are suppression and 
enhancement of normal diffusion. Our results are explained by a simple theoretical approximation. 
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To understand diffusion in noisy maps, that is, in time-discrete dynamical systems where the deterministic equations 
of motion are perturbed by noise, figures as a prominent problem in recent literature. The most simple example of such 
models are one-dimensional chaotic maps on the line. In seminal contributions by Geisel and Nierwetberg and by 
Reimann et al. scaling laws have been derived for the diffusion coefficient yielding suppression and enhancement 
I . of diffusion with respect to variation of the noise strength. Related results have been obtained in Refs. However, 
all these results apply only to the onset of diffusion where the scaling laws are reminiscent of a dynamical phase 
transition, and not much appears to be known far away from this transition point. In such more general situations, 
■ only perturbations by a nonzero average bias have been studied |^]. Related models are deterministic Langevin 
\ equations, in which the interplay between deterministic and stochastic chaos has been analyzed however, without 
. focusing on diffusion coefficients. Non-diffusive noisy maps have furthermore been investigated by refinements of cycle 
expansion methods Q . 

In this work we study the transition scenario from deterministic to stochastic diffusion in the most simple type of 
chaotic dynamical systems, which are piecewise linear maps on the line. Particularly, we are searching for signatures 
of deterministic and stochastic dynamics in the diffusion coefficient as a function of the strength of time-dependent 
stochastic noise. In this aspect our work appears to be related to the recent dispute on a possible distinction between 
^<r^ \ chaotic and stochastic diffusion in experiments 1^ , where some of the theoretical models studied are very similar to 
ps) . the one introduced below. 

' We define our system as follows: The unperturbed map is given by the equation of motion 
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Xn+l = Ma{Xn) , (1) 



, where a S M is a control parameter and x„ is the position of a point particle at discrete time n. Ma{x) is continued 
^ ' periodically beyond the interval [—1/2, 1/2) onto the real line by a lift of degree one, Ma{x -I- 1) = Ma{x) + 1. We 
assume that Ma{x) is anti-symmetric with respect to x = 0, Ma{x) = —Ma{—x). The map we study as an example is 
^ ■ defined by Ma(x) = ax, where the uniform slope a serves as a control parameter. The Lyapunov exponent of this map 
\ is given by A = In a implyingthat for a > 1 the dynamics is chaotic. We now apply two types of annealed disorder to 
this map, (i) noisy slopes 1^,^: we add the random variable Aa„, n £ N, to all slopes a making them time-dependent 
^\ in form of 

Ma+Aa„ {x) = (a + Aa„)a; , (2) 

or (ii) noisy shifts H-l: we add the random variable A6„, n G N, as a time-dependent uniform bias yielding 

MaMx) = ax + A6„ . (3) 

In both cases we assume that the random variable A„ G {Aa„,A6„} is independent and identically distributed 
according to a distribution Xd{^n)i where d G {da, db} is again a control parameter. In the following we will consider 
two different types of such distributions, namely random variables distributed uniformly over an interval of size [—d, d] 

xd(A„) = ^e(d + A„)e(d-A„) , (4) 

and dichotomous or (5-distributed random variables 12,01 , 
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XdiAn) = ^{S{d ~ An) + Sid + A„)) . (5) 

Since |A„| < d, we denote d as the perturbation strength. As an example, we sketch in Fig. ^ our model for noisy 
slopes. We now define the diffusion coefficient as 

D{a,d) = lim -!-(< xl >p„ - < a;„ >^ ) , (6) 

n— too Zn 

with 

< a;^ >po= J dx J rf(Ao)d(Ai) . . . d(A„„i)po(a;)x(Ao)x(Ai) • ■ • x(A„-i)a;^ , (7) 

where poix) denotes the initial distribution of an ensemble of moving particles, Xo = x, k N, and Aj,j £ 
{1, . . . ,n — 1}, is the random variable. Note that in computer simulations it suffices to generate a single series of 
random variables instead of evaluating all the integrals in Eq. (|^) . To obtain better numerical convergence for noisy 
shifts the current squared in Eq. was subtracted at any time ste p w hile Eqs. (^),(||) imply that the long-time 
average over the random variable A„ does not yield any bias. In Refs. [ |To| , pj] | it was shown that the unperturbed map 
Eq. (|l|) exhibits normal diffusion if a > 2, and the same was found recently by adding a bias b [|l2|. Correspondingly, 
for the types of perturbations defined above diffusion should always be normal if (a — da) > 2, as was confirmed in 
simulations. Hence, the central question is what happens to the parameter-dependent diffusion coefficient D[a,d) 
under variation of the two control parameters a and d in case of the above two types of noise. 

For da = it was shown that the unperturbed diffusion coefficient D(a, 0) is a fractal function of the slope a as a 
control parameter |^,^, as is depicted again in Fig. ^. Included are results from computer simulations for uniformly 
distributed noisy slopes at different values of the perturbation strength da ||l3|. As expected, the fractal structure 
gradually smoothes out by increasing da. Qualitatively the same result is obtained by applying noisy shifts [ p^ . 
Fig. |l| may be compared to the corresponding result for quenched slopes Fig. 1 in Ref. Apart from numerical 
uncertainties, there are clear differences in the critical behavior close to the onset of diffusion. However, for small 
enough perturbation strength and large enough a the results look qualitatively similar indicating that in this limit 
quenched and annealed diffusion may be treated on the same footing. 

More evidence for this statement is obtained from a trivial approximation for the perturbed diffusion coeffi- 
cient, which we motivate starting from dichotomous noisy slopes. Naive reasoning suggests that, at arbitrary 
fixed parameters a and da, the perturbed diffusion coefficient D{a, da) can be approximated by simply averag- 
ing over the unperturbed diffusion coefficient D{a,0) at respective values of the slopes a — da and a + da yielding 
Da,pp{a,da) = {D{a — da, 0) -I- D{a + da, 0))/2. This heuristic argument can be straightforwardly extended to the 
random distribution Eq. as well as to any other type of uncorrelated noise yielding the generalized expression 

Dapp(p, d) = y d(A) Xdi^)Dip + A, 0) . (8) 

Here p is a vector of control parameters such as p = {a, 6} in case of the map above, d is the corresponding vector of 
perturbation strengths, and A is the vector of perturbations such as A = {Aa, Ab} for noisy shifts and slopes. Further 
generalizations of this equation, for example, to arbitrary moments as defined in Eq. are straightforward. Applying 
this formula to the case of quenched slopes discussed in Ref. j|] reproduces the diffusion coefficient approximation Eq. 
(6) therein, which was obtained in the limit of small perturbation strength. The corresponding approximations for 
uniform noisy slopes are depicted in Fig. |l| as lines. They show that even for the rather large perturbation strength 
da = 1 the agreement between theory and simulations is excellent. This confirms that, in the limit described above, 
quenched and annealed disorder generating normal diffusion can indeed approximately be treated in the same way. 

Let us now look at the diffusion coefficient for a given value of a as a function of da. Fig. ^ shows that approxi- 
mately at odd and even integer slopes the fractal diffusion coefficient D{a,0) exhibits a local maximum or minimum, 
respectively. Since Eq. (H) represents an average over the unperturbed solution in a local environment [a — da, a + da] it 
predicts local suppression and enhancement of diffusion at odd and even integer slopes, respectively, under variation 
of the perturbation strength da. This has already been conjectured in Ref. and has been verified in Ref. j|] 
for quenched slopes. We first check this hypothesis for noisy slopes around the local maximum of D(a,0) at a 7 
distributed according to Eqs. (@),(||). Figs.|| (a), (b) depict again results obtained from computer simulations in com- 
parison to Eq. (|^). As predicted, in both cases there is suppression of diffusion for small enough da. For dichotomous 
noise the perturbed diffusion coefficient increases on a coarse scale by exhibiting multiple, fractal-like suppression 
and enhancement on finer scales. For uniform perturbations there is a pronounced crossover from suppression to 
enhancement on a coarse scale, by again exhibiting oscillations on a fine scale. In both cases the agreement between 



2 



simple theory and simulations is excellent for small enough da, whereas clear systematic deviations particularly in 
case of dichotomous noise are visible for larger da. Note that if a — Aa^ < 2 particles are getting trapped within a box 
at a respective time step n, and that for a — Aa„ < 1 the map is non-chaotic. In the first case simulations and simple 
reasoning suggest that the perturbed map still exhibits normal diffusion. However, as soon as a — da < 1 numerical 
results indicate that there is no normal diffusion anymore |l4[ . This appears to be due to the contracting behavior of 
the non-chaotic map resulting in localization of particles. The oscillatory behavior of the diffusion coefficient in Fig. 
^ (a) just below this transition point is not yet understood. 

Employing Eq. (||) we now analyze noisy shifts. The unperturbed two-parameter diffusion coefficient D{a,b,0) has 
been calculated numerically exactly for the map under consideration in Ref. . Results for the perturbed diffusion 
coefficient D{a, db) = D{a, 0, db) are presented in Fig. ^ (a) for dichotomous noise and (b) for uniform perturbations, 
both starting from D(a,0) at a = 6. In both cases the perturbed diffusion coefffcient exhibits strong enhancement of 
diffusion for small enough perturbation strength due to the fact that the unperturbed diffusion coefficient at a = 6 
is approximately identical with a local maximum in the (a, 6) parameter plane p^ . For dichotmous perturbations 
it suffices to show results for < d& < 0.5 only. Translation and reflection symmetry of the map imply that this 
function is mirrored in the interval from 0.5 < dfe < 1, and that the full sequence in < d6 < 1 is periodically repeated 
for higher values of db. As in the corresponding case of noisy slopes, the perturbed diffusion coefficient increases on 
a coarse scale by exhibiting multiple fractal-like suppression and enhancement on a fine scale. In case of uniform 
perturbations there is a pronounced crossover to an approximately constant diffusion coefficient for larger db. 

Before calculating the stochastic limit of the diffusion coefficient we provide a simple analytical justification for the 
heuristic approximation Eq. (H). For sake of simplicity, we demonstrate it only for noisy slopes, A„ = Aa„. Noisy 
shifts as well as quenched disorder can be treated along the same lines [Q. Let us start from the definition of the 
diffusion coefficient Eq. where < x„ >= 0. Let Aa„ be uniformly distributed in [—da, da], Auq = Aa. In case of 
da all random variables are bounded by Aa„ = Aa -I- e, —2da < e < 2da. We now put this expression into the 
perturbed equation of motion Eqs. (|l|),(|^), as contained in Eq. (^, which we write as Xn+i,a+Aa„ = Ma+Aa^ixn)- As 
a first step we now take the limit e ^ resulting in the expression for the mean square displacement 



<xl> = I dxj d{Aa)d{Aai) . . .d{Aan-i)po{x)x{Aa)x{Aai) . . .x{Aan-i)xl^^^^^^_^ 

dx J diAa)poix)xiAa)xl,+^^ (e 0) . (9) 
As a second step we exchange the time limit contained in Eq. with the integration over d{Aa) yielding 

Dapp{a, da) lim 
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d{Aa)xda{Aa) lim / dxpo{x 
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n,a+Aa 
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= J d{Aa)xda{Aa)D{a + Aa, 0) , (10) 
where we have used that the unperturbed diffusion coefficient was defined as 

D{a,0)= lim dxpo{x)xl (11) 

n — ^oc / 

We have thus verified our previous approximation Eq. (||) for noisy slopes in the limit of small perturbation strength. 
A similar derivation can be carried out for noisy shifts arriving again at Eq. (^ in case of very small perturbation 
strength. For quenched shifts it is known that a normal diffusion coefficient does not exist [ p^ , thus any approximation 
by Eq. (H) must fail. Indeed, it turns out that in this case taking the limit e — > fundamentally changes the properties 
of the dynamical system and is thus no valid operation ]T^ . 

Finally, we calculate the parameter-dependent stochastic diffusion coefficient related to the map with noisy slopes. 
Starting from the definition Eq. (11) the complete loss of memory in the unperturbed map is modeled by 16 



(i) replacing the distance Xn & particle travels by n times the distance a particle travels at any single time step, 
nAx = n{Ma{x) — x), and (ii) neglecting any memory effects in the probability density on the unit interval by 
assuming po{x) — 1. Then Eq. ( |ll| ) yields 

(a - I'l^ 

DUa) = . (12) 
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As was shown in Refs. ]ll| , p^ , this equation correctly describes the asymptotic parameter dependence of the determin- 
istic diffusion coefficient for a — > cx) thus explaining the increase of D{a, 0) in Fig. |l] on a coarse scale. On this basis, 
the corresponding result for noisy slopes is easily calculated by using Eq. (|l^) as the functional form for D(a + Aa, 0) 
in the approximation Eq. (^) reading 

Drw{o., da) — Drw{a,0) + ^a'^ /c , (13) 

where c = 24 for dichotomous noise Eq. (|^) and c ~ 72 for uniform noise Eq. (|). Eq. (|l|) thus confirms the common 
sense expectation that noise should typically enhance diffusion and represents the stochastic limit of the diffusion 
coefficient. This equation is depicted in Fig. || (a), (b) in form of dashed lines. In case of dichotomous noise the 
correlations are apparently large enough such that even for large perturbation strength da there is no transition to 
the stochastic limit, whereas in case of uniform noisy slopes the diffusion coefficient approaches the stochastic solution 
asymptotically in da thus verifying the existence of a transition from deterministic to stochastic diffusion. That such 
a distinct transition behavior exists in these models was already conjectured in Ref. Analogous calculations 

for noisy shifts yield Eq. ( [l^ for all values of db reflecting the fact that for large enough a the stochastic diffusion 
coefficient should not depend on the bias. This result is shown in Fig. |^ (b) and again confirms an asymptotic approach 
of the diffusion coefficient to the stochastic limit under variation of db. Based on the known result of the existence of 
a fractal diffusion coefficient for the unperturbed D{a, b, 0) we conjecture that the typical transition scenario in this 
type of systems consists of (multiple) suppression and enhancement of diffusion. We finally note that Eqs. (|l2|),(p^) 
are closely related to the approximation outlined in Ref. , and to the simple heuristic argument given by Reimann 
1^ by which he explains the suppression of deterministic diffusion by noise in the climbing sine map near a crisis; 
more details will be discussed elsewhere p4[ |. 

We conclude with a few remarks: (1) It would be interesting to study the problem of noisy maps with non-zero 
average bias along the same lines. Ref. [ p^ shows that the unperturbed map does not exhibit linear response for 
6 — > 0, thus we conjecture that adding noise generates a transition to Ohm's law. (2) In the recent Ref. a 
rather general mechanism of noise suppression by noise has been reported. Whether there is a more detailed relation 
between the argument outlined in this reference and the phenomena discussed here appears to be an open question. 
(3) Our approach may be useful to investigate the impact of noise on the diffusion coefficient in more complex time- 
continuous systems as well. In particular, we are thinking of models such as the standard map, particle billiards, or 
inertia ratchets, where irregular transport coefficients have already been reported and studied under the impact of 
noise p^ ]. However, these analyses were not performed from the point of view of suppression and enhancement of 
diffusion, or by looking for transitions to the stochastic limit. 

The author thanks N. Korabel and S. Dcnisov for interesting discussions on noisy maps. He is also grateful to G. 
Radons, H. van Beijeren, J.R. Dorfman, and T. Tel for helpful remarks. 
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FIG. 1. Diffusion coefficient D{a,da) for the piecewise linear map shown in the figure. The slope a is perturbed by uniform 
noise of maximum strength da as defined in Eq. (^). The bold black line depicts numerically exact results for the unperturbed 
diffusion coefficient at da — 0. Computer simulation results for da ^ are marked with symbols, the corresponding lines 
are obtained from the approximation Eq. (^). The parameter values are: da = 0.1 (circles), da = 0.4 (squares), da = 1.0 
(diamonds). 
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FIG. 2. Diffusion coefficient D{a, da) as a function of the perturbation strength da at slope a = 7 for noisy slopes distributed 
according to: (a) dichotomous noise Eq. (|^, (b) uniform noise Eq. (^. The circles represent results from computer simulations, 
the bold lines are obtained from the approximation Eq. (^), the dashed lines represent the stochastic limit for the diffusion 
coefficient Eq. (p^). 




FIG. 3. Diffusion coefficient D{a, db) as a function of the perturbation strength db at slope a = 6 for noisy shifts distributed 
according to: (a) dichotomous noise Eq. (b) uniform noise Eq. (^. The circles represent results from computer simulations, 
the bold lines are obtained from the approximation Eq. (^), the dashed line represents the stochastic limit for the diffusion 
coefficient Eq. (p^). 
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